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Abstract. We present a new approach to derive the observed 
properties of synchrotron outbursts in relativistic jets. The idea 
is to use the very well sampled submillimetre-to-radio long- 
term light curves of 3C 273 to extract the spectral and temporal 
evolution of a typical outburst. The method consists in a de- 
composition of these light curves into a series of twelve self- 
similar flares. With a model-independent parameterization, we 
find that the obtained outburst's evolution is in good qualitative 
agreement with the expectations of shock models in relativistic 
jets. We then derive, by a second approach, the relevant param- 
eters of three-stage shock models. 

We observe for the first time that the optically thin spectral 
index is steeper during the initial rising phase of the evolution 
than during the final declining phase as expected by the shock 



model of Marscher & Gear (1985). We obtain that this index 



flattens from a = — 1.1 {Si, oc ^") to a = —0.5, in good agree- 
ment with what is expected from a power law electron energy 
distribution of the form N{E) ex E^"^-^. The observed flatten- 
ing gives support to the idea that radiative (synchrotron and/or 
Compton) losses are the dominant cooling process of the elec- 
trons during the initial phase of the shock evolution. 

Two other results give us confidence in our decomposi- 
tion: 1) the outbursts that we identify do well correspond to 
the VLBI components observed in the jet and 2) there is strong 
evidence that high-frequency peaking outbursts evolve faster 
than low-frequency peaking outbursts. We propose that this last 
correlation is related to the distance from the core of the jet at 
which the shock forms. 

Key words: galaxies; active - galaxies: jets - quasars: individ- 
ual: 3C 273 - radio continuum: galaxies 



1. Introduction 

Submillimetre-to-radio light curves of blazars show evidence 
of prominent structures, or flares, apparently propagating from 
high to low frequencies. A decisive step in the understanding 
of these flares was done by Marscher & Gear (1985, hereafter 



MG85). They studied the strong 1983 outburst of 3C 273 by 
constructing at two epochs a quasi-simultaneous millimetre- 
to-infrared spectrum after subtracting a quiescent emission as- 
sumed to vary on a much longer time scale. They successfully 
fitted these two flaring spectra with self-absorbed synchrotron 
emission and showed that their temporal evolution can be un- 
derstood as being due to a shock wave propagating down a rel- 
ativistic jet. They identified three stages of the evolution of the 
shock according to the dominant cooling process of the elec- 
trons: 1) the Compton scattering loss phase, 2) the synchrotron 
radiation loss phase and 3) the adiabatic expansion loss phase. 
Another shock model was developed by Hughes et al. 



(1985) simultaneously to that of MG85. Their piston-driven 
shock model reproduces well the lower frequency flux and po- 
larization observations of outbursts in BL Lacertae, but fails 
to describe the observed behaviour in the millimetre domain. 
A generalization of the three-stage shock model of MG85 was 



presented by Valtaoja et al. (1992). Their model, based on ob 



servations, describes qualitatively the three stages of the MG85 
model without going into the details of the physics of the shock. 



Finally, Qian et al. (1996) proposed a burst-injection model to 



study the spectral evolution of superluminal radio knots. Their 
theoretical calculation is able to reproduce well the ob served 
spectral evolution of the C4 knot in 3C 345 (Qian |l996| ). 

To constrain these shock models, we need to extract the 
properties of the outbursts from the observations. This step is 
difficult both at high and at low frequencies. At high frequen- 
cies because of the brevity of the outbursts that last only a few 
days to months, thus requiring a very well sampled set of ob- 
servations in the not easily accessible submillimetre spectral 
range. At radio frequencies because they very often overlap due 
to their longer duration, making it difficult to isolate them. 

The best observational constraints for the model of MG85 



were obtained by Litc hfield et al. (|1995|) for the blazar 3C 279 
and by Stevens et al. (|I995[ Il99al998|) for PKS 0420-014, 
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3C 345 and 3C 273, respectively. All these studies are based on 
isolated outbursts. The method used consists in constructing si- 
multaneous multi-frequency spectra for as many epochs as pos- 
sible after the subtraction of a quiescent spectrum assumed to 
be constant with time. The subtraction of a quiescent spectrum 
is convenient and seems to give good results, but has only weak 
physical justification. In 3C 273, there was a period of nearly 
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constant flux at millimetre frequencies lasting just more than 
one year in 1989-1990, which was interpreted as its quiescent 
state (Robson et al. 1993). At radio frequencies, however, no 



similar constant flux period was ever observed and there is no 
evidence that such a state exists at a significant level above the 
contribution of the jet's hot spot 3C 273A (see Fig. 2 of Tiirler 
et al. 



1999ai 



The different approach presented here to derive the ob- 
served properties of the outbursts has the advantage to not rely 
on the assumption of a quiescent emission. The idea is to de- 
compose a set of light curves covering a large time span into 
a series of flares. To our knowledge, the first attempt of such 



a decomposition was made by Legg (1984), who fitted a ten 



years radio light curve of 3C 120 with twelve self-similar out- 
bursts. Recently, Valtaoja et al. ( 1999[ ) decomposed the 22 GHz 
and 37 GHz radio light curves of many active galactic nuclei 
into several exponentially rising and decaying outbursts. What 
is new in our approach is that we fit the same outbursts simulta- 
neously to twelve light curves covering more than two decades 
of frequency from the submillimetre to the radio domain. This 
adds a new dimension to the decomposition: the evolution of a 
flare is now a function of both time and frequency. The aim is 
to obtain both the spectral and temporal properties of a typical 
flare, from which individual flares differ only by a few param- 
eters. 

We use the light curves of 3C 273, the best observed quasar, 
to have as many observational constraints as possible. The flar- 
ing behaviour of 3C 273 was alrea dy the subject of seve ral pre - 
vious studies (e.g. Robson et al. 1993; Stevens et al. 1998). 



Stevens et al. (1998) obtain results for the first stage of the 
strong 1995 flare in very good agreement with the predictions 
of the MG85 shock model. The new approach presented here is 
however more powerful to constrain the two following stages 
of the evolution. 

We describe below two different approaches. In Sect, g we 
model the light curve of each outburst by an analytic function 
that can smoothly evolve with frequency, whereas in Sect. H 
we directly model a self-absorbed synchrotron spectrum that 
evolves with time. The first approach is easier to implement, 
since it allows us to begin the decomposition with a single light 
curve before adding the others progressively. The second ap- 
proach is more physical and gives better constraints to shock 
models. Our results are discussed in Sect. &and summarized in 
Sect. |. 

Throughout this paper the frequency v is as measured in 
the observer's frame and "log" refers to the decimal logarithm 
"logj^o "■ The convention for the spectral index a is 5iy oc v 



+ Q 



2. Observational material 

This study is based on the fight curves of the multi-wavelength 



database of 3C 273 presented by Tiirler et al. (1999c). The 



twelve light curves we use are the five radio light curves: 
5 GHz, 8.0 GHz, 15 GHz, 22 GHz and 37 GHz and the seven 
millimetre/sub millimetre (mm/sub mm) light curves: 3.3 mm, 
2.0 mm, 1.3 mm, 1.1mm, 0.8 mm, 0.45 mm and 0.35 mm. At 



low frequency (5 to 15 GHz), we consider only the measure- 
ments of the University of Michigan Radio Astronomy Ob- 
servatory (UMRAO). The observations at 22 GHz and 37 GHz 
are mainly from the Metsahovi Radio Observatory in Fin- 
land. The mm/submm observations are from various sources 
including the James Clerk Maxwell Telescope (JCMT), the 
Swedish-ESO Submillimetre Telescope (SEST) and the "Insti- 
tut de Radio- Astronomic Millimetrique" (IRAM). 

We analyse the observations from 1979.0 to 1996.6, except 
at low frequency where we extend the analysis up to: 1997.2 
(15 GHz), 1997.5 (8.0GHz) and 1998.0 (5 GHz), in order to 
include the decay of the 1995 flare. In the mm/submm range, 
we average repeated observations made within 3 days to avoid 
oversampling of the light curves at some epochs. This leaves 
us a total of 4352 observational points to constrain our fits. To 
observations without known uncertainties, we assign the aver- 
age uncertainty of the other observations at the same frequency. 
The light curves are treated as if all their observations were 
made exactly at the same frequency, i.e. small differences of 
the observing frequency from one measurement to the other 
are not taken into account. This simplification should not much 
affect the results, since the spectrum is rather flat (a <L — 0.5) 
in the con sidered submillimetre-to-radio domain (e.g. Tiirler et 
al. 



1999a) 



3. The light-curve approach 

We describe here an approach in which we minimize the num- 
ber of model-dependent constraints. The light curve of each 
outburst at a given frequency is described by a simple analyt- 
ical function. The choice of this function is purely empirical 
and does not rely on any physical model. The evolution with 
frequency of the outburst's light curve is left as free as possi- 
ble. This model has therefore many free parameters, which can 
adapt to a wide range of different situations. 

3.1. Number of outbursts 

One crucial parameter of the decomposition is the number of 
outbursts. Pushed by the wish to reproduce the small features 
seen in the light curves, one is tempted to add always more 



outbursts to the fit. In Tiirler et al. (1999b), we published the 



results of a decomposition into nineteen outbursts using an ap- 
proach which is similar to that described below. Here, we try 
to minimize as much as possible the number of outbursts to 
better constrain their spectral and temporal evolution. We end 
with twelve flares, which are absolutely necessary to describe 
the main features of the light curves. 

The aim of the decomposition is not to reproduce the de- 
tailed structure of the light curves, but to derive the main char- 
acteristics of the outbursts. As a consequence, the x^ that we 
shall obtain will be statistically completely unacceptable and 
will have no meaning in terms of the probability that the model 
corresponds to what is observed. We will however refer to the 
obtained values of the reduced y^ (cf. Sect. 3J), because it is 
the usual way to express the quality of a fit. 
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Fig. 1. Model light curve of an outburst defined by Eqs. ([l|) 
and (g]), which starts at time to = and peaks at an amplitude 
of A{i') — 1 Jy. The three different line types show the effect of 
varying p{i') and (/)(zy) 



3.2. Parameterization 

At a given frequency v, we model the light curve S^{t) of a 
single outburst of amplitude A{v), starting at time t — to and 
peaking att = to+ iriso(;^) by 



S.{t) 



Aiu) 



t-tp 

tiisc{v) 



pH' 



ifto <t <to + triso(i^) and by 



S^{t) = A(i^)exp 



t-to-triscjiy) 



0MN 



(1) 



(2) 



if t > to + tiiscii^)- The exponents p(i^) and 0(z^) define the 
shape of the light curve at frequency i' and tfan(i^) is the e- 
folding decay time of the flare at frequency i/. Different time 
profiles of an outburst defined by Eqs. (lib and (0) are shown in 

Fig.[i|. 

Rather than constraining the outburst parameters iA{i'), 
triscii'), tia\i{i^), p{v) and 4){v)) at each of the twelve light 
curve's frequencies, we describe their logarithm by a cu- 
bic spline which we parameterize at only four frequencies 
spaced by 0.75 dex and covering the 3 -600 GHz range (see 
Fig. H). This reduces the number of free parameters by a factor 
three, while keeping the parameterization completely model- 
independent. We thus need a total of 5 x 4 parameters to fully 
characterize the spectral and temporal evolution of an out- 
burst, i.e. a surface in the three dimensional (5, v, t)-space (cf. 

Fig. a. 



We impose that all individual outbursts are self-similar, in 
the sense that they all have the same evolution pattern, i.e. the 
same shape of the surface in the (S*, v, t)-space. What we al- 
low to change from one outburst to the other is the normaliza- 
tion in flux 5, frequency v and time t, which changes, respec- 
tively, the amplitude of the outburst (strong or weak), the fre- 
quency at which the emission peaks (high- or low-frequency 
peaking) and the time scale of the evolution (long-lived or 
short-lived). A change in normalization corresponds to a shift 
of the position of the outburst's characteristic surface in the 
(logs', log i^, log t) -space. To define this position, we take the 
point of maximum flux as an arbitrary reference point on the 
surface. On average among all individual outbursts, this point 
is located at ( (log S) , (log v) , (log t) ) and this average normal- 
ization defines what we call the typical outburst of 3C 273. 
We denote by AlogS, Alogj/ and Alogi the logarithmic 
shifts of this point with respect to the average position, i.e. 
A log A: = logfc - (log A;), 'ik = S,v,t. These 12 x 3 loga- 
rithmic shifts plus the 12 different start times to of the flares 
give a total of 48 parameters used to define the specificity of all 
outbursts. 

The superimposed decays of the outbursts that started be- 
fore 1979 are simply modelled by an hypothetical event of am- 
plitude Ao{v) at time t = to + ti.isc(i^) = 1979.0 and de- 
caying with the c-folding time tfaii(i^) of the typical outburst 
at frequency v. The variation of the amplitude Ao{v) with 
frequency is modelled by a cubic spline as for the five other 
variables, but parameterized at four slightly lower frequencies 
(log(z^/GHz) == 0.5, 1.0, 1.5 and 2.0), due to the fact that 
Ao{v) is only well constrained for the radio light curves. Fi- 
nally, we assume a constant contribution to the light curves due 
to the quiescent emission of the jet's hot spot 3C 273A. This 
emission is modelled with a power law spectrum as given in 
Turler et al. (|l999a|). 



To summarize, this first parameterization uses a total of 72 
(20 + 48 + 4) parameters to adjust the 4352 observational points 
in the twelve light curves. The great number of free parameters 
still leaves more than four thousand degrees of freedom (d.o.f.) 
to the fit. The simultaneous fitting of the twelve light curves 
is performed by many iterative fits of small subsets of the 72 
parameters. 

3.3. Results 

Fig. H shows three representative light curves among the twelve 
fitted simultaneously with the outbursts parameterized as de- 
scribed in Sect. 3.2. The major features of the light curves are 
reproduced by the model with only about one outburst every 
1.5 year starting simultaneously at all frequencies. The overall 
fit has a reduced x^ value of x^cd — X^ /'^■'^■^■~^^-^- The main 
discrepancy between the model and the observations arises dur- 
ing 1984-1985, when the very different light curve features in 
the millimetre and radio domains cannot be correctly described 
by the 1983.4 flare alone. 

The obtained evolution of the parameters with frequency 
for the typical outburst is shown in Fig. g. The amplitude 
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Fig. 2. a-c. Three of twelve light curves fitted by the cumu- 
lative sum (solid line) of twelve outbursts (dotted lines) pa- 
rameterized with the light-curve approach described in Sect. ^. 
a The 1.1mm (- 280 GHz) light curve (x?od = 30.5); b the 
37 GHz hght curve (x?ed = 21-4); c the 8.0 GHz light curve 
(Xred = 9-6)- Each outburst starts simultaneously at all frequen- 
cies at the epoch Iq shown in panel a. The dashed line is the 
contribution of the jet's hot spot 3C 273A 



A{iy) of the light curve has a maximum at ~ 45 GHz. Both 
the rise time irise(i^) and the c-folding decay time ifaii(i^) in- 
crease monotonically with wavelength. If we extrapolate the 
cubic spline to low frequency, it is striking to see that both 
tiisc{i^) and <fau(j^) tend to very high values of the order of 
10 years at 1 GHz, while the amplitude of the outburst would 
still be significant (^(1 GHz) « 1 Jy). Due to the lack of sub- 
millimetre observations before 1981, the amplitude Ao{i') is 
not constrained at frequencies above ^ 300 GHz. The increase 
of AQ{h') at these frequencies - due to the spline - is proba- 
bly not real, but does not affect the fit because the correspond- 
ing decay time is short (ifaii( 1000 GHz) « 1 year). The two 
exponents p{iy) and (p{iy) which describe the shape of the out- 
burst's light curve are both higher at radio frequencies than in 
the mm/submm domain. As a consequence, the light curves at 
higher frequencies have a steeper rise just after the start of the 
outburst and a steeper decay just after the peak (see Figs. |l] and 
0a). 

The five parameters A(i^), triso(j^), itan(i^), p{i^) and (^{i') 
define the typical outburst that can be represented in three di- 
mensions in the (log S, log ly, log t)-space as shown in Fig. 0d. 
The three other panels of Fig. show the three Cartesian pro- 
jections of this surface. The frequency and time axes cover the 
same logarithmical range of 4 dex, so that the dotted diagonal 
in Fig. ^ corresponds to v xt^^. At least at low frequencies. 
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Fig. 3. Evolution with frequency of the six parameters: A(iy_ 
Aoiiy), t,-isc{t^), ttaii(j^), p{t^) and <j){i^), defined in Sect. 3.2. 



All these functions are a cubic spline passing through the four 
points with frequencies fixed at log (I'/GHz) — 0.5, 1.0, 1.5 
and 2.0 for ^o(t^) and at log (i^/GHz) ~ 0.5, 1.25, 2.0 and 
2.75 for the five other parameters 



both the maximum of the spectra and of the light curves follow 
quite well this diagonal. The outburst's evolution is thus amaz- 
ingly symmetric in Figs, ^a and ^. The maximum amplitude 
of the typical outburst is of ^ 16 Jy and is reached after ^^ 7.5 
months at a frequency of ~ 45 GHz. The frequency i^ni of the 
spectrum's maximum is steadily decreasing with time (Fig. ^). 
The corresponding flux density Sm is first increasing with de- 
creasing frequency i/m according to 5m oc i^m"'^, whereas it 
decreases as 5m oc ly^^'^ during the final decline of the outburst 
(Fig. ^). This behaviour corresponds qualitatively to what is 
expected by shock models (e.g. MG85). 

At frequencies above the spectral turnover (v :^ ;-',„), the 
spectral index a is first of ^ —0.5 and steepens very slightly 
to ^ —0.7 at the maximum development of the outburst. The 
somewhat chaotic behaviour during the final declining phase 
- due to the abrupt change in the parameter 0(i/) - does not 
enable us to define a reasonable spectral index during this last 
stage. At frequencies below the spectral turnover {i/ -^ i/m), the 
spectral index a is smoothly steepening with time from '-^ 1.8 
to ~ 2.5. This is what is expected from a synchrotron source 
that starts inhomogeneous and progressively becomes homo- 
geneous (e.g. Marscher 1977 1. 

The twelve individual outbursts have different amplitudes 
ranging from 5 Jy up to 32 Jy for the 1983.0 flare studied by 
MG85. The corresponding dispersion a of the amplitude shifts 
A log 5 is (T = 0.23, which is slightly smaller than the disper- 
sion of the frequency shifts Alogz^ (a — 0.27) and the time 



M. Tiirler et al.: Modelling the submillimetre-to-radio flaring behaviour of 3C 273 



Time [ year 
0.01 0.1 1 



Frequency [ GHz ] 
1 10 100 1000 




0.01 0.1 1 10 

Time [ year ] 



^og(u) 



Fig. 4. a-d. Logarithmic spectral and 
temporal evolution of the typical out- 
burst obtained by the light-curve ap- 
proach described in Sect. |[ Panel d 
shows the three-dimensional represen- 
tation in the (log 5*, log ly, log t)-space. 
The other panels show the three Carte- 
sian projections: a light curves at dif- 
ferent frequencies spaced by 0.2 dex; 
b spectra at different times spaced by 
0.2 dex; c contour plot in the frequency 
versus time plan. The thick solid and 
dashed lines show the evolution of the 
maximum of the spectra and of the light 
curves, respectively. They follow quite 
well the dotted diagonal in panel c cor- 
responding to 1/ oc t~^. The star dot 
shows the point of maximum develop- 
ment of the typical outburst. The open 
circles show for each outburst the posi- 
tion that would have this maximum ac- 
cording to the shifts A log S, A log i^ 
and Alogi. The arrows in panel b 
show the frequency distribution of the 
twelve light curves 



shifts A log t (a — 0.29). The amplitude shifts A log S are ob- 
viously not correlated with either A log t/ or A log t (Fig. ^ 
and b). This is confirmed by a Spearman rank-order test (Bev- 
ington 1969), which yields that the observed correlations could 
occur by chance with a probability of more than 60 %. On the 
contrary, the shifts A log i/ and A log t align well along the 
V oc t~^ line (Fig. ^) and the Spearman's test probability of 
< 0.01 % confirms that this anti-correlation is very significant. 

4. The three-stage approach 

In the light-curve approach described above, we model analyt- 
ically the light curve of an outburst at different frequencies and 
show that the resulting typical flare is qualitatively in agree- 
ment with what is expected by shock models in relativistic jets. 
It is thus of interest to derive from the data the parameters that 
are relevant to those models. 

The sho ck mo del of MG85 and its generalization by Val- 
taoja et al. ( 1992 ) describe the evolution of the shock by three 
distinct stages: 1) a rising phase, 2) a peaking phase and 3) a 
declining phase[|. The three-stage approach presented below is 
similar to that of Valtaoja et al. (1992 ), in the sense that its aim 
is simply to qualitatively describe the observations. It contains 



We use here the terminology introduced by Qian et al. ( 1996 ), 
because it is purely descriptive and free of any interpretation regarding 
the physical origin of these stages. 



however more parameters in order to include those which are 
relevant to test the physical model of MG85. 



The remarks of Sect. 3.1 concerning the number of out- 
bursts and the quoted values of the reduced x^ ^PPly equally 
here. 



4.7. Parameterization 

The self-absorbed synchrotron spectrum emitted by electrons 
with a power law energy distribution of the form N{E) oc E^^ 
can be expr essed - by generalizi ng the homogeneous case (e.g. 



Pacholczyk 



197C; Stevens etal. 



1995) -as 



Si, — S\ 



"""^ 1 - exp(-(i^/i/i)"""»~"""='') 



1-c- 



(3) 



where (iy/iy]^)"thin-athick is equal to the optical depth t^, at fre- 
quency ly. Si and i/i are respectively the flux density and the 
frequency corresponding to an optical depth of Ti, = 1 . At high 
frequency (ly ^ vi) the medium is optically thin (t^, ^ 1) and 
the spectrum follows a power law of index athin^^ (s — l)/2, 
whereas at low frequency (i/ <|; i/i) it is optically thick (ti, ^ 1) 
and the spectral index is athick- In the case of a homogeneous 
source, athick = +5/2. 

The maximum 5m = 5'^(fm) of the spectrum S,y is reached 
at the turnover frequency f^ corresponding to an optical depth 
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of Tm — (i^m/'^i)""''" "thick ^^^ jg obtained by differentiating 
Eq. %. 

^ = ^ exp(r„0-l=fl-^iii^)r„,. (4) 

Av \ athick / 



By developing the exponential of Eq. ^ 
obtain a good approximate: r,n = \\\j 



to the third order, we 
We 



can now rewrite Eq. 

r,n and ^m by 



Sathin _ 2 
3 athick 

according to the turnover values v^ 



Du — On 



V 



1 - exp (-Tm {v/v^Y 



1 



(5) 



The evolution with time of the self-absorbed synchrotron 
spectrum of Eq. (Q) is assumed to follow three distinct stages: 
1) the rising phase for t — io < ^r ; 2) the peaking phase for 
ir < t — to < ip and 3) the declining phase for t — to > tp. 
The subscripts "r" and "p" refer to the end of the rising phase 
and the end of the peaking phase, respectively. We assume that 
during each stage i (i = 1, 2, 3) both the turnover frequency 
Vmit) and the turnover flux 5m (i) evolve with time as a power 
law, but with exponents that differ during the three stages: 



i(i)cxt''' and S'mli) oc F' 



Sn. oc ly^/f"' 



(6) 



We thus need ten parameters: t,-,tp, fm(ir), Sm{tr), Pi, (^2, Pa, 
7i, 72 and 73, to describe the evolution of the spectral turnover 
in the three dimensional {S, v, i)-space. 

The model of MG85 predicts that both the optically thin 
Qfthin and thick athick spectral indices should be flatter during 
the declining phase than during the rising and peaking phases 
(see Fig. 3 of Marscher et al. 1992). To test whether the spec- 
trum is actually changing from the rising phase to the declining 
phase, we allow the two spectral indices athin and atinck to 
have different values during these two stages. The transition 
during the intermediate peaking phase from the values in the 
rising phase (Q;thin(ir) and athick(ir)) to the values in the de- 
clining phase (Q;thin(ip) and athick(ip)) is assumed to be linear 
with the logarithm of time log (i). This adds the four parame- 
ters athin(*r), Q;thin(ip), Q;thick(ir) and athick(ip) to the model, 
having thus a total of fourteen parameters to fully define the 
evolution of a typical flare in the {S, v, t)-space instead of the 



twenty parameters used in the first approach (Sect. 3.2) 



The specificity of each outburst is modelled with a total 
of 12 X 4 parameters exactly as described in Sect. 3.2 for the 
light-curve approach. We do not model again the superimposed 
decays of the outbursts that started before 1979, but simply 
use the same exponential decay as obtained by the first ap- 
proach (Sect. p!^ ). The constant contribution of the jet's hot 
spot 3C 273A is also considered here. The total number of pa- 
rameters in this second parameterization is a bit less than for 
the first one: 62 (12 x 4 + 14) instead of 72. 

4.2. Results 

To allow a better comparison with the results of the first ap- 
proach (Sect. 3^), we show in Fig. Q|the same light curves as 



X) 20 




Fig. 5. a-c. Same as Fig. |[ but with the outbursts parameter- 
ized according to the three-stage approach described in Sect. ^ 
a The 1.1mm (^ 280 GHz) light curve (x^^^ ^ 26.8); b the 
37 GHz light curve (x?od = 20.4); c the 8.0 GHz light curve 



ixt 



10.9) 



Table 1. Values of the parameters defined in Sect. \^_A_ 



responding to the evolution of the typical outburst shown in 
Fig. H The two first columns display the fourteen best fit pa- 
rameters, whereas other related parameters are shown in the 
last column 



Param. 


Value 


Param. 


Value 


Param. 


Value 


ir 


0. 14 year 


ip 


1.63 year 






fm(ir) 


120 GHz 






l'm(ip) 


13.8 GHz 


Sra{t,) 


15.3 Jy 






Oin(ipj 


15.9 Jy 


athin(ir) 


-1.09 


athin(ip) 


-0.48 






athick(ir) 


+ 1.55 


athick(ip) 


+ 1.74 






/3i 


-0.51 


71 


+0.51 


71 //3i 


-0.99 


A 


-0.88 


72 


+0.02 


72//32 


-0.02 


A 


-1.19 


73 


-1.36 


73//33 


+ 1.14 



in Fig. g. The reduced x^ of the overall fit is now of x^ed = 17.8. 
The higher frequency light curves are relatively better de- 
scribed here than with the first approach (compare Figs, g and 
||). The start times ig of the outbursts are very similar to those 
obtained by the first approach, except for the fourth flare which 
is now starting much later at tg = 1984. 1 instead of 1983.4. This 
later ip seems to be in better agreement with the observations, 
but the behaviour of 3C 273 during 1984-1985 is still poorly 
described. 

The obtained values of the parameters are given in Table |l[ 
They correspond to the spectral and temporal evolution of the 
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Fig. 6. a-d. Logarithmic spectral and 
temporal evolution of a typical flare in 
3C 273 obtained by the three-stage ap- 
proach described in Sect. H. The range 
covered by the axis is the same as in 
Fig. H. All lines and points are defined 
as described in Fig. H. The maximum 
of the spectrum follows the thick soUd 
line. The three different slopes of this 
line correspond from left to right in 
panel a to 71, 72 and 73; in panel b to 
73 //?3, 72 //32 and 71 //3i; and in panel 
c to /3i, (32 and (3^. The values of these 
slopes are given in Table [l|. Notice in 
panel b the change of the optically thin 
spectral index athin from —1.1 during 
the rising phase to —0.5 during the de- 
clining phase. This change is responsi- 
ble for the strange shapes of the high 
frequency light curves in panel a 



typical outburst shown in Fig. g. If the tracks followed by the 
maximum of the spectra and of the light curves are similar 
to those obtained by the first approach (Fig. ^, the spectral 
evolution of the outburst derived here is quite different. We 
obtain that the spectral turnover flux Sj^ increases during the 
first 50 days (t^ — 0.14 year) with decreasing turnover fre- 
quency z^in as ^in oc i'^^'^ ■ The subsequent very flat peak- 
ing phase is found to be relatively long, since it lasts 1 .5 year 
and spans nearly one order of magnitude in frequency from 
120 GHz to 13.8 GHz. The final declining phase is quite abrupt 
with a relation between S„i and t^m of S^ oc t'^^'^. The opti- 
cally thin spectral index athin is found to be clearly steeper in 
the rising phase than in the declining phase. It is flattening by 
Aathin = +0.6 during the peaking phase from athin (^r ) = — 1 • 1 
to athin(ip) = —0.5. The optically thick spectral index athick 
is found to be more constant with a slight tendency to steepen 
with time. It has a mean value of athick = +1.65 and is steep- 
ening by Aathick — +0.2 during the peaking phase. 

For each outburst we obtain logarithmic shifts in amplitude 
A log S, frequency A log i' and time A log t, which are simi- 
lar to those obtained by the first approach (Sect. BJ). The dis- 
persions (T of A log S, A log ly and A log t are 0.20, 0.34 and 
0.27, respectively. A possible correlation of A log S with either 
A log 1/ or A log t is again not significant: the Spearman's test 
probability that stronger correlations could occur by chance is 
> 40 %. On the contrary, the strong correlation observed be- 



tween A log u and A log t is most probably real (Spearman's 
test probabihty < 10"'^). 



5. Discussion 

The two approaches presented above give comparable results, 
but differ concerning the existence of a nearly constantpeak- 
ing phase and the shapes of the spectra (compare Figs. ^ and 
^). The origin of these differences can be understood by com- 
paring the light curve profiles obtained by the two approaches, 
which are shown in Fig. 0. It is clear that the first approach al- 
lowing only a rising phase and a declining phase cannot mimic 
the three-stage profiles of Fig. |7]b resulting from the second 
approach (Sect. 4.1). On the other hand, only the light-curve 
approach is able to produce a round peaking phase as seen in 
Fig. 0a. In a forthcoming paper (Tiirler et al. in preparation), 
we will present the results of an hybrid approach, which in- 
corporates the advantages of both approaches in order to better 
define the properties of the typical outburst. 



5.7. Do the outbursts correspond to VLBI components? 

The decomposition of the light curves into distinct outbursts 
was motivated by the observation with very long baseline inter- 
ferometry (VLBI) of distinct components in the jet structure of 
3C 273. Since the detection of a new VLBI component (Krich- 
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Time [ year 



Fig. 7. a and b. Light curves at different frequencies for the 
typical outburst obtained with the light-curve approach (a) and 
the three-stage approach (b). The dotted light curves are spaced 
by 0. 1 dex in frequency and the six solid light curves are at fre- 
quencies log (i^/GHz) = 3.0, 2.5, . . . , 0.5, in order of increas- 
ing time scales 



baum et al. 1990 1 associated with the strong optical/infrared 



flare of 1988 in 3C 273 (Courvoisier et al. |1988[ ), there is good 
evidence that outbursts are related to the ejection of new VLBI 
knots. To test whether all outbursts are actually associated with 
superluminal components, we compare in Table the start time 
ip of an outburst - as obtained by the three-stage approach - 
with the ejection time "tg (knot)" of a new VLBI knot as given 



by Abraham et al. (1996 1 and Zensus et al. (1990). For each 
of the eight first outbursts, we can identify one or two possibly 
associated VLBI components. 

To test further this relationship, we compare the flux densi- 
ties "i^obs (knot)" of the VLBI components observed at epoch 
t = 199L15 and at a frequency of 10.7 GHz (Abraham et al. 
1996() with the flux densities Fcxp(t = 1991.15, i/= 10.7 GHz) 



expected at the same epoch and the same frequency according 
to the outburst parameters derived here. Table || shows that for 
the five first outbursts there is always one of the possibly as- 
sociated knots (indicated by an arrow), which has the expected 
flux. For the three remaining outbursts and especially for the 
1988.1 flare, the relation between Fobs (knot) and F^xp is not 
obvious. At this epoch however, the possibly associated com- 
ponents are still strongly blended by the core emission (com- 
ponent "D") or might even still be part of the unresolved coreg. 
The total flux Foxp = 20.9 Jy expected by the 1986.3, 1988.1 
and 1990.3 outbursts is indeed equal to the observed total flux 

^ In our model the core emission is entirely due to a superimposition 
of outbursts unresolved by the VLBI. 



Table 2. Relation between VLBI components and the eight 
first outbursts as obtained with the three-stage approach 
(Sect. H). The parameters in this table are defined in Sect. 5.1 



to 



Knot 



to (knot) 



Fc: 



Fobs (knot) 



1979.6 


C5 


1978.6 ±0.04 


0.4 Jy 


0.7±0.4Jy ^ 




C6 


1980.0 ±0.04 




1.9±0.5Jy 


1980.9 


C6 


1980.0 ±0.04 


1.9 Jy 


1.9±0.5Jy ^ 




X 


< 1985.6 




<0.5Jy 


1982.4 


C7 


1982.2 ±0.4 


0.2 Jy 


<0.5Jy ^ 


1983.1 


C7a 


1983.1 ±0.00 


0.9 Jy 


<0.5Jy 




C7b 


1983.6 ±0.09 




1.2±0.3Jy ^ 


1984.1 


C7b 


1983.6 ±0.09 


3.8 Jy 


1.2±0.3Jy 




C8 


1984.7 ±0.10 




4.4±0.5Jy ^ 


1986.3 


Cx 


< 1988.2 


1.0 Jy 


<0.5Jy 


1988.1 


C9 


1988.4 ±0.17 


11.2Jy 


1.8±0.4Jy 


1990.3 


CIO 
D 


< 1990.2 


8.7 Jy 


6.4±0.5Jy 
13.2 ± 0.3 Jy 



Fobs = 21.4 ± 0.7 Jy of the C9, CIO and D components. These 
results strongly suggest that there is a close relation between 
the outbursts and the VLBI knots and hence that our decompo- 
sition describes a real physical aspect of the jet. 

5.2. How can we understand the peculiarities of individual 
outbursts? 

The relation found between the outbursts and the VLBI knots 
(Sect. 5.1) has established that our decomposition is not 
purely mathematical, but does correspond to a physical real- 
ity. There should therefore be a physical origin to the clear 
anti-coiTelation found between the frequency shifts Alogi^ 
and the time shifts Alogi of the individual outbursts. The 
observed frequency shifts A log i' confirm that 3C 273 emits 
both low- and high-frequency peaking outbursts (Lainela et 
al. 1992). The relation between Alogj/ and Alogt clearly 



shows that high-frequency peaking flares evolve faster than 
low-frequency peaking outbursts. The alignment of the shifts 
along the v oc t^^ line (Figs. Eb and He) further suggests the 
relation A log t/ = — A log t. 

The origin of this relation could be due to a change A log V 
of the Doppler factor V = 7~^(1 — /3cos9)^^, which de- 
pends on the flow speed /3 = v/c, the Lorentz factor 7 = 
(1 — /3^)^^/^ and the angle to the line of sight 6. Observed 
quantities (unprimed) are related to emitted quantit ies (p rimed) 
as (e.g. Hughes & Miller |l 99 ll Pearson & Zensus |l987| ): 



SH 



t^v-^t' 

= P3 s'{v') 



A log J/ = A log P + A log v' (7) 

Alogt = -AlogP + Alogt' (8) 
AlogS-^SAlogP + AlogS" (9) 



If we assume that in the jet frame all outbursts are alike (i.e. 
A log k' = Q,\/ k = S, ly, t), the observed relation A log i/ = 
—A log t can be interpreted as a change A log V of the Doppler 
factor from one outburst to the other. In this case, however. 
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there should also be correlations between A log S and both 
A log I' and A log t, which are not observed. 

Alternatively, we can consider that the Doppler factor does 
not change (AlogP — 0) and that the observed relation be- 
tween A log ly and A log t is intrinsic and independent of pos- 
sible flux variations A log S. Such a correlation might be re- 
lated to the distance from the core at w hich th e shock forms 
(Lainela et al. |1992|). Indeed, Blandford (|1990|) shows that for 



a simple conical jet with constant speed v the frequency of 
maximum emission j/m is inversely proportional to the dis- 
tance down the jetr — vt (i^m ex r~^), while the corresponding 
flux density S'm is constant. Since the speed v is constant, the 
turnover frequency j/m is then also inversely proportional to 
time (i/jji c>c t~^), as observed. If a shock forms in such an un- 
derlying jet at a distance tq from the core, both the frequency 
range of the emission and the time scale of the evolution will 
depend on the distance rg, as illustrated in Fig. 0. We therefore 
propose that short-lived and high-frequency peaking flares are 
actually inner outbursts, whereas long-lived and low-frequency 
peaking flares are outer outbursts. 

This interpretation is supported by the existence of short- 
lived VLBI components which are only seen close to the core. 
In our decomposition, the two most short-lived and the most 
high-frequency peaking outbursts are the two successive flares 
of 1982.4 and 1983.1. Their start times correspond well to the 
period from 1981 to 1983 during which only short-lived VLBI 



components were formed (Abraham et al. 1996). If our inter- 
pretation is right, the shifts A log i/ and A log t that we obtain 
suggest that the 1982.4 flare would have formed about two 
times closer to the core than the 1983.1 flare and four times 
closer than the typical outburst. 

5.3. What are the constraints for shock models? 

According to the shock model of MG85, the optically thin 
spectral index athin should be steeper during the two first 
stages of the outburst evolution than the usual value of athin = 
— (s — l)/2 (Sect. 4.1). A steeper index arises due to the fact 



that the thickness x of the emitting region behind the shock 
front is proportional to the cooling time icooi of the electrons 
suffering radiative (Compton and/or synchrotron) losses. Dur- 
ing the rising and peaking phases, radiative losses are domi- 
nant and therefore the thickness x is frequency dependent as 
X oc icooi oc i^~^/^, which leads to a steeper optically thin spec- 
tral index of atinn = — s/2. Until now, the expected flattening 
of the spectral index by Aathin = +0.5 from the rising and 
peaking phases to the declining phase was never observed and 
furthermore the optically thin spectral index athin observed at 
the beginning of the outburst was often found to be already too 
flat (athin > — 1/2) to allow t he exp ected subsequent flattening 
(Valtaoja et al. |1988| Lainela |1994 . 

The present result that the optically thin spectral index is 
flattening with time by Aathin = +0.6 is in good agreement 
with the change of Aatinn — +0.5 expected by the shock model 
of MG85. The observed flattening of the spectrum is contrary 
to the steepening with time expected as a result of radiative en- 




ar. 



Fig. 8. Schematic representation of a conical jet with con- 
stant speed V in which the frequency of maximum emission v^ 
(greyscale) is inversely proportional to the distance down the 
jet r = tit (i^m oc r^^ oc t^^). Let us assume that the peak- 
ing phase of an outburst corresponds to an octave of radius 
(r -^ 2 r). This phase has a duration and a frequency range 
which clearly depend on the distance tq at which the shock 
forms. Inner outbursts peaking from r*/2 to r* are four times 
shorter than outer outbursts peaking from 2r^, to 4r*, while 
their emission is maximum at four times higher frequencies 



ergy losses by the electrons. The observed behaviour can how- 
ever also be understood as a change of slope with frequency 
rather than with time and thus it could conceivably be due to a 
spectral break that steepens the optically thin spectral index by 
a factor of 0.5 at higher frequencies. Such a break is expected in 
the case of continuous injection or reacceleration of electrons 
suffering radiative losses (Kardashev 1962) and is observed 
in several hot spots including 3C 273A (Meisenheimer et al. 
1989). Whatever the interpretation, the flatter index, athin(ip). 



is the relevant index to determine that the electron energy in- 
dex s (N{E) oc ^-'') is s = 1 - 2 athin(ip) = +2.0. This value 
corresponds to the a verage value observed in several hot spots 
(Meisenheimer et al. 1989) and is in agreement with the values 
expected if the electrons are acceler ated b y a Fermi mechanism 
in a relativistic shock (e.g. Longair 1994). 

The long flat peaking phase observed in 3C 273 contrasts 
with t he co mplete absence of this stage in 3C 345 (Stevens 
et al. 1996). This difference is surprising, because the out- 



burst's evolution is otherwise very similar in these two objects 
with nearly the same indices for the rising and the declining 
phases: 7i//3i = -0.99 in 3C 273 and -0.86 in 3C 345 and 

73//33:=+1.14in 3C 273 and +0.98 in 3C 345 (S'^oci/^i'^'^')- 
A value of 73 //33 ~ +1 was also found in several other sources 



by Valtaoja et al. (1988). This decrease of the turnover flux 



with decreasing frequency is steeper than expected by the sim- 
plest model of MG85; i.e. with a conical adiabatic jet having 
a constant Doppler factor T). With s == 2 and a magnetic field 
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B oriented perpendicular to the jet axis, their model predicts 
I3/P3 = +0.45. This discrepancy between the observations and 
the shock model of MG85 was akeady pointed out by Stevens 



et al. ( 1996 ). We refer the reader to their discussion of two more 



general cases of the MG85 model: 1) a straight non-adiabatic 
jet and 2) a curved adiabatic jet. With the observed values of the 
indices f]^ and 73, these authors could determine the two free 
parameters of the model. In our case, with the constraints of all 
six indices /3i and 7^ (z = 1, 2, 3), we could not find a good 
agreement with either of the two models mentioned above. In 
a forthcoming paper (Tiirler et al. in preparation), we will fur- 
ther discuss this point and explore whether a non-conical non- 
adiabatic curved jet can well describe the observations. 



6. Summary and conclusion 

By using most available submillimetre-to-radio observations of 
3C 273, we have been able to extract the properties of the spec- 
tral and temporal evolution of a typical outburst. The new ap- 
proach we defined consists in decomposing the light curves into 
several self-similar outbursts. The main results of our decom- 
position are the followings: 

- It is possible to understand the very different shapes of the 
submillimetre-to-radio light curves of 3C 273 with only 
about one outburst every 1.5 year starting simultaneously 
at all frequencies. 

- There is no need to invoke any underlying quiescent emis- 
sion apart from the weak contribution of the jet's hot spot 
3C 273A. 

- The outbursts that we identify do well correspond to the 
observed VLBl components in the jet. 

- There is good evidence that short-lived and high-frequency 
peaking flares are emitted closer to the core of the jet than 
long-lived and low-frequency peaking outbursts. 

- The spectral and temporal evolution of the outbursts is 
found to be in good qualitative agreement with the evolu- 
tion expected by shock models in relativistic jets. 

- We observe a flattening of the optically thin spectral index 
from the rising to the declining phase of the shock evolu- 
tion, which supports the idea proposed by MG85 that ra- 
diative (synchrotron and/or Compton) losses are the main 
cooling process of the electrons during the initial phase of 
the outburst. 

We are aware that our decomposition is far from describing 
the detailed structure of the light curves and that the jet emis- 
sion is much more complicated than this work tries to show. 
Nevertheless, the results suggest that the outbursts we iden- 
tified are closely related to the VLBl knots, and hence that 
they describe a physical aspect of the jet. The new approach 
presented here is a powerful tool to derive the observed prop- 
erties of millimetre and radio outbursts. It allows comparison 
between shock models and the observations and we are con- 
fident that such decompositions are able to further constrain 
present and future shock models. Finally, we would like to 



stress the importance of long-term multi-wavelength monitor- 
ing campaigns, which turn out to be essential towards a better 
understanding of the physics involved in relativistic jets. 
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